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We study the effects of spatially localized breakdown of time reversal symmetry on the surface of a 
topological insulator (TI) due to proximity to a helical spin density wave (HSDW). The HSDW acts 
like an externally applied one-dimensional periodic(magnetic) potential for the spins on the surface 
of the TI, rendering the Dirac cone on the TI surface highly anisotropic. The decrease of group 
velocity along the direction x of the applied spin potential is twice as much as that perpendicular 
to x. At the Brillouin zone boundaries (BZB) it also gives rise to new semi-Dirac points which 
have linear dispersion along x but quadratic dispersion perpendicular to x. The group velocity of 
electrons at these new semi-Dirac points is also shown to be highly anisotropic. Experiments using 
TI systems on multiferroic substrates should realize our predictions. We further discuss the effects 
of other forms of spin density wave on the surface transport property of topological insulator. 



PACS numbers: 71.10.Pm, 73.20.-r 

I. INTRODUCTION 

In the last few years there has been a growing interest 
in topological insulators (TI), which are materials that 
are insulators in the bulk but conduct two-dimcnsionally 
on the surfaces^— . The non-trivial topology of the wave- 
function of such TI has been predicted to show metal- 
lic surface conductivity that is topologically protected 
against weak disorder and interaction effects. One of the 
mechanisms by which the surface modes can be disrupted 
is by breaking the time reversal invariance upon applying 
a magnetic field. There have been several studies on the 
effects of a magnetic field cither applied directly perpen- 
dicular to the surface of a Tl£ or by the proximity effect 
of a fcrromagnct in a hcterostructure^^. On the other 
hand, the effects of an antiferromagnet or a spin density 
wave on the surface states of the TI is a topic which is 
much less well-explored. We carry out such an analysis 
of the effect of a spin density wave on TI transport in 
this work. 

Although an antiferromagnet, in the proximity to a TI, 
does not affect the global time reversal symmetry (TRS) 
in the latter, the staggered nature of the spins in the an- 
tiferromagnet breaks the TRS locally. In this work, we 
are interested in how this local TRS breaking affects the 
surface states of the TI. However, the local TRS breaking 
occurs at a length scale which is of the order of the lattice 
spacing in the antiferromagnet (for e.g., 8.85A for MnO). 
Therefore the effects of local TRS breaking occur at very 
large momenta, which makes it difficult to observe them 
in experiments. In order to redress this difficulty, we pro- 
pose to study the effects of local TRS breaking by replac- 
ing the antiferromagnet with a multiferroic material 
like orthorhombic RM11O3 (R being a rare earth element 
like Tb or Dy) or the family of Fei-^Co^Si with cubic 



but noncentrosymmetric structure, which shows a helical 
spin density wave (HSDW) order with relatively long pe- 
riods (> 30nm)ii. The HSDW will couple to the Dirac 
fermions on the surface of the topological insulator due 
to the proximity effect by the above two kinds of ma- 
terials directly deposited on the top. The situation is 
similar to that of an externally applied charge potential 
on graphenei^. However, the contrast between the two 
situations is that the 3D topological insulator has spin 
polarized 2D Dirac fermions on its surface^ 3 -, whereas 
the Dirac fermions in graphene carry pseudo-spin. Nev- 
ertheless, in analogy with the charge potential applied to 
graphene, a spin potential can be used for TI to manip- 
ulate the surface Dirac fermions. 

It may be appropriate to ask about the motivation of 
our theoretical work where we are proposing, using con- 
crete theoretical calculations, that experiments be car- 
ried out on a multifcrroic-TI sandwich structure for the 
observation of the topologically protected surface trans- 
port in the TI surface states induced by the HSDW asso- 
ciated with the multiferroic layer. This approach sounds 
somewhat indirect, and indeed it is although the effect 
of the HSDW-induced local TRS-breaking on TI surface 
transport properties is intrinsically interesting in its own 
right as we establish in this work. 

Our main motivation for this work, however, arises 
from the fact that so far there has been little direct ex- 
perimental signature of the topologically-protected sur- 
face transport in existing TI materials in spite of a great 
deal of research activity-^— aimed precisely toward the 
direct and unambiguous observation of 2D surface TI 
transport. All the convincing experimental evidence for 
the existence of the topologically-protected surface 2D 
TI states is currently based on the verification of the pro- 
posed band structure through STM/STS or ARPES type 
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spectroscopic measurements. The problem in the direct 
observation of the surface 2D TI transport is that the cur- 
rently existing TI materials are not good bulk insulators, 
and the bulk conduction is always much stronger than the 
surface conduction, making it impossible to observe the 
surface 2D states unambiguously in transport measure- 
ments. This problem of substantial bulk transport arises 
from the fact that all existing TI materials, instead of be- 
ing a bulk band insulator as band structure calculations 
predict, turn out to be intrinsically doped in the bulk due 
to defects and vacancies, giving rise to a large bulk con- 
duction channel which competes directly with the surface 
2D transport. For example, the two most recent (and also 
most compelling) transport measurement o 14 ^ 15 see very 
small putative surface 2D magneto-resistance oscillations 
with a temperature dependent bulk resistivity which docs 
not behave like a standard band insulator at all. In par- 
ticular, the bulk resistivity in these sample s 14 ' 15 is of the 
order of O-cm or less whereas an insulator typically has a 
bulk resistivity which is 10-12 orders of magnitude larger. 
Thus, the TI systems, even in these most compelling 
measurements, show dominant bulk conduction with less 
than 1% of the net conduction being inferred (through 
the indirect fitting of the magneto-resistance data using 
many free parameters and ad hoc multichannel conduc- 
tion models) at best to be arising from the 2D surface 
states. Another recent experiment^ concludes that no 
2D surface conduction can be discerned at all in the TI 
transport data because of the dominant bulk conduction. 
Another recent experiment^ goes to the extreme of us- 
ing a pulsed external magnetic field as high as 55T in 
order to investigate the surface 2D transport, again em- 
phasizing small observed features in the high-field resis- 
tivity as the possible manifestation of the expected 2D 
TI surface transport. Even the two very recent transport 
measurement o 14 ' 15 purportedly claiming the manifesta- 
tion of 2D surface TI transport can only observe small 
2D features in the derivative spectrum of the resistiv- 
ity with respect to the applied magnetic field. This is a 
most unsatisfactory state of affairs in sharp contrast to 
the transport properties of well-established 2D quantum 
systems, e.g. graphene, where 2D transport behavior— 
without any bulk conduction problem whatsoever mani- 
fests itself in every possible transport measurement in a 
decisive manner, and one does not have to look for small 
features in the derivative spectrum. The current situa- 
tion in TI physics is thus extremely problematic with all 
spectroscopic measurements providing reasonable verifi- 
cation of the expected spin-resolved 2D Dirac cone band 
structure on the TI surface whereas the fact that the 
bulk is not an insulator is making it essentially impossi- 
ble to see the expected 2D protected surface transport. 
We point out in this context that the theoretical details 
of how the surface 2D TI transport would behave, had 
it not been contaminated by the unintentional bulk con- 
duction problem plaguing the existing TI materials, are 
reasonably well-known in the literature^. 

Of course, further materials development in the exist- 



ing TI systems leading to the effective suppression of the 
unintentional bulk doping or the discovery of completely 
new classes of TI materials where the bulk is a true insu- 
lator could solve this problem instantaneously, but until 
that happens, any idea which points toward the obser- 
vation of surface TI transport should be welcome. It is 
somewhat of an embarrassment that in spite of the huge 
activity in TI research, a true Topological Insulator does 
not yet exist since the existing systems are not true bulk 
insulators due to the invariable presence of unintentional 
bulk doping. Our proposal in the current work should be 
seen in this light. We provide a method to directly isolate 
surface transport features in TI systems in proximity to 
a helical spin density wave as produced, for example, by 
a multiferroic material. Since bulk conduction should be 
relatively immune to the presence of the HSDW, whereas 
the surface TI conduction should be affected qualitatively 
as we show in this work, we believe that our predictions 
could help provide the unambiguous observation of the 
surface 2D TI transport properties. In addition, the cou- 
pling between the HSDW and the 2D TI states leads to 
nontrivial novel physics (e.g. the semi-Dirac points to be 
discussed below) which is intrinsically interesting in its 
own right. 

We find that the surface states of the TI are sig- 
nificantly modified in the presence of the HSDW. The 
HSDW acts like an externally applied superlattice poten- 
tial on the TI surface resulting in striking anisotropy of 
the Dirac cones and group velocity of the surface states. 
In particular, the group velocity along the direction of 
the superlattice is monotonically suppressed as a func- 
tion of the lattice potential strength and its period. For 
the type of HSDW considered in Sec. HH which we call 
proper HSDW^ 4 - (see below), we find that novel semi- 
Dirac points, whose low-energy characteristics are inter- 
mediate between Dirac (masslcss) and zero-gap (massive) 
scmiconductor o 25 ' 26 , show up on the BZB. In Sec. IIII1 
we discuss how the surface transport properties of the 
topological insulator changes when the forms of the spin 
density wave is changed. Section ITVl contains a summary 
and conclusions. 



II. THE PROPER HELICAL SPIN DENSITY 
WAVE 

A. Dirac cone on the BZ center 

We begin by writing down the effective Hamiltonian for 
the low-energy quasiparticlcs on the surface of a topolog- 
ical insulator as, 

H (k) = hv F (k x a x + kydy) (1) 

where vf — 6.2 x 10 5 m/s is the Fermi velocity of Dirac 
fermions in Bi2Se;j2i and Oi are the Pauli matrices. The 
Hamiltonian in Eq. ([1]) is characterized by an energy 
spectrum e s .k = shvpk where s — ±1 is the band index, 



and cigenstatcs given by 
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where 6*k is the angle of vector k with respect to the k x 
direction. 

In the presence of the HSDW on top of a TI, depicted 
schematically in Fig. [IJ the Hamiltonian (QJ is modified 28 
to: H = Ho(k) +U(x), where the potential U(x) for the 
HSDW can be written as: 
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U{X) = UyOy COS( — X 
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U z a z sin(—x), (3) 



where L is the spatial period of the potential, and U V:Z are 
the amplitudes of the HSDW—. For the case U z < U y , we 
find no gapless states on the BZB. On the other hand, for 
Uz > U y , two Dirac points centered at (j^,±K y ) emerge 
on the BZB, which will be discussed in Sec. IIIII For 
the symmetric HSDW such that U z = U y = Uo, which 
we shall call proper HSDW, we find a single semi-Dirac 

point at (k x ,k y ) = (± — ,0) (Fig. [3^,). In this section, we 

shall focus on the case of a proper HSDW on TI. 

In our calculation, we use the value of induced ex- 
change field due to the magnetic proximity effect being 
5 ~ 50 me V 29 ' 30 , which has been taken to be a reason- 
able valued. The exchange- induced potential U(x) from 
the HSDW creates a periodic potential on the surface of 
the TI. The electronic eigenstates of such a Hamiltonian 
can be written using Bloch's theorem as 
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x (m £ integer) is the reciprocal lattice 

vector, ip n k{r) and c(n, k + G) are 2-spinor functions and 
n is the band-index. Here k x is only limited to be in the 

first Brillouin Zone (FBZ), i.e. k x e [--, -]. The band 
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eigenstates ip n k(f) and the corresponding eigen-energies 
E n (k) are obtained as eigenvalues and eigenvectors of the 
Bloch equation 

[hv F ((k x + G)a x + kycj y ) - E nk ]c(n,k + G) 

2ti\ 



(U y a v + iU z a z )c{n, k + G- — )/2 

Lj 

{UyOy - iU z a z )c(n, k + G + ^)/2 



(5) 



While the above Bloch equation may be solved by nu- 
merical diagonalization to obtain numerical eigenvalues 
and eigenvectors, more insight can be obtained into the 
solutions around high symmetry points in the FBZ at 
k = and at the BZB by using perturbation theory. 
The numerical results of diagonalizing Eq. [5] are shown in 
Fig. 2, Fig. 3 and Fig. 4. The total Hamiltonian H Q + U 
of the TI+HSDW system is invariant under a composite 
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FIG. 1: (Color online) (a) Schematic diagram of the system. 
The HSDW is on top of the topological insulator with the 
periodicity L. (b) Energy of charge carrier dispersion versus 
the wavevector k with k measured near the BZ center. Dashed 
(green) and solid (blue) lines show the dispersion along the 
x direction and along the y direction, respectively with Uo = 
50meV and L — 90nm (£7(0, 0) was shifted to zero). 



symmetry S = 8T where O is the time-reversal opera- 
tor and T is the operator corresponding to translation by 
L/2. Since, the operator S, similar to the time-reversal 
symmetry operator O, is both anti- unitary and satisfies 
S 2 = —1, the proximity to the HSDW docs not open a 
gap at fc = 0. In other words, the HSDW preserves the 
global time-reversal symmetry, it does not open up a gap 
near fc = on the TI surface. 

In the limit of a perturbatively weak coupling to the 
HSDW ( fanncfl ^ I); t ne leading order contribution of 

0) dispersion can 
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the HSDW to the low-momentum (fc 
be characterized by the renormalization of the group ve- 
locity and effective mass at fc = 0. The group velocity of 
states without the HSDW is isotropic around k = with 
constant magnitude vp. In presence of the HSDW, the 
renormalized group velocity of quasiparticles parallel to 
the wavevector k [vf, = v(k) • k] around the Dirac point 
(| A; I <C 1) obtained within second order perturbation the- 
ory can be written as 



V k ^ V F 

vf 



I d(E+, k -e +tk ) 



flVF 



dk 



UqL 2 
8ir 2 h 2 v 2 F 



(3 + cos 20 k ) 
(6) 

From Eq. ((6j), it is clear that the renormalized group ve- 
locity is anisotropic around the Dirac point and decays 
monotonically, both with the amplitude of periodic po- 
tential Uo, and with the spatial period L of the potential. 
Fig- HI shows the result from a full numerical calculation. 
We note that there is good agreement between the trends 
from Eq. © and the full calculation (shown in Fig. 
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FIG. 2: (Color online) (a) The group velocity Vf, (measured 
from the Dirac point) of charge carriers on the surface of the 
TI in units of the Fermi velocity {vf) versus the angle (9k 
with L — 90nm. Dotdashed red, solid green and dashed blue 
lines correspond to Uo being 10 meV, 30 meV and 50 meV, 
respectively, (b) The inverse effective mass me/m*^ versus 
Uo for the first band above the Dirac point on BZ center. 
Dotdashed red, solid green and dashed blue lines are results 
for L being 30nm, 60nm and 90nm, respectively. 



for weak potential strength. 

The interaction with the HSDW introduces at k ~ 
a finite curvature along the ^-direction to the previously 
linear in k dispersion of the Dirac fermions. The effective 
mass tensor at k = from the curvature of the energy 
band around the Dirac point is given by 
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is the bare electron mass and Go = — . From 

L 

Eq. (JT)) the effective mass m e /m* x grows as U(j, which 
describes the small Uq behavior of the numerically deter- 
mined effective mass shown in Fig. ^jp. 



B. Dirac Cone near BZ boundary 

While the periodic potential from the HSDW in- 
troduces renormalizations of the velocity and effective 
masses at the center of the Dirac cone at fc = 0, the 
effect of the periodic potential is limited to being pertur- 
bative since there is no degeneracy in the initial Dirac 
spectrum near k = 0. However, such degeneracies do oc- 
cur at the edges of the FBZ, which we have referred to 
as the BZB, and non-perturbative effects of the potential 
may be generated by coherent back scattering with mo- 
menta G. Unlike the case of conventional crystals, where 
such back-scattering opens up gaps, one finds the emer- 
gence of new Dirac cones at the BZB in the presence of 
interaction with the HSDW (Fig. [3Ji). 

When the wavevector k is on the first BZB ((k x , k y ) = 
(ir/L, k y )), the two states |s, k) and |s, k — (2n/L, 0)) are 
degenerate before the periodic potential is applied. In the 
presence of an applied potential, the largest contribution 
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FIG. 3: (Color online) (a) Semi-Dirac charge carrier disper- 
sion versus the wavevector k with k measured from (k x , k y ) = 

TV 

(y,0) on the BZB. Dashed (green) and solid (blue) lines 

Li 

show the linear dispersion along the x direction and the 
quadratic dispersion along the y direction, respectively with 
U = 50meV, L = 90nm and Go = 27r/90nm _1 (E(n/L,0) 
was shifted to zero). (b).The energy gap AE between the 
first and the second band at the BZB versus k y for charge 
carriers above the Dirac point with Uo = 50meV and Go = 
27r/90nm~ 1 . Dotdashed (red) , solid(green) and dashed (blue) 
lines correspond to L being 30 nm, 60 nm and 90 nm, respec- 
tively. 



to the energy eigenvalues at the edges of the BZ comes 
from these two degenerate states. For the clockwise helix, 
the backscattering amplitude leads to an energy gap on 
the BZB that is given by, 



AE(k y ) = Uo 1 - 



y/1 + (kyL/ny 



(8) 



where s — ±1 represents the positive or negative bands 
in the Dirac cone. In the lower band of the Dirac cone, 
i.e. s = — 1, the above gap is maximum at k y = and 
decreases monotonically with \k y \. The positive band is 
more interesting with a gap that vanishes at k y = 0. 
From the full numerical calculation in Fig. we can 
see that the energy gap of the positive band on the BZB 
increases monotonically then decreases with k y . 

If we change the sign of both U y and U z , the physics 
will not change since the total Hamiltonian is invariant 
under this transformation T, where T is the operator 
corresponding to translation by L/2 in the x direction. 
The chirality of the helix depends on the relative sign of 
the U y and U z in the HSDW potential. The effects of 
the clockwise helix on the upper band of the Dirac cone 
is the same as the anticlockwise helix on the lower band 
of the Dirac cone. For the anticlockwise helix, the minus 
sign in Eq. [5] should be changed to a plus sign, which 
leads to a vanishing gap at k y = for the lower band. 

As with the Dirac cone at k = 0, an understanding of 
the low-energy dispersion and the corresponding trans- 
port properties associated with the Dirac cone at the 
BZB is obtained by calculating the group velocity and 
effective mass. Second order perturbation theory leads to 
a strongly anisotropic group-velocity vt measured from 
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FIG. 4: (Color online) (a) Anisotropic Vf,/vF measured from 
the semi-Dirac point on the BZB versus the angle #k with 
L = 90nm. Dotdashed (red), solid (green) and dashed (blue) 
lines correspond to Uo being 10 meV, 30 meV and 50 meV, 
respectively, (b) The inverse effective mass m e /m* y versus 

TV 

L around the point (k x ,k y ) — (y,0) for the lower band of 

the semi-Dirac point for three Uo- Note the negative effective 
mass (corresponding to Fig. [3] (a)). 



the gap-less point at {k x ,k y ) = ( — ,0) which is written 
as 



hiipGo cos 6*k 
JU$ + h 2 v\G 2 



(9) 



For comparison, the full numerical calculation is shown 
in FiglU The effective mass tensor at the new gapless 
point for the clockwise helix, calculated within second 
order perturbation theory similar to Eq. is given by 



0. 



k=(ir/L,0) 



2(Hv F G ± U )/Gl a,p = y 

(10) 

where the ± corresponds to the upper or lower band 
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around the semi-Dirac point at (k x ,k y ) = ( — ,0). Eq. 

IH1 and QJJ together show the highly anisotropic dispersion 
of the semi-Dirac point on BZB: a dispersion linear along 
the periodic direction x but quadratic along the perpen- 
dicular direction y. For the anticlockwise helix, the ef- 
fective mass tensor given in Eq. [TUJ is still valid except 
now the gap for the lower band vanishes at k y = 0. Fig. 
EH shows the angular dependence of the group velocity 
^ic from a full numerical calculation. The velocity v y is 
seen to have zero component regardless of the strength 
of the applied potential. The dependence of the renor- 
malized group velocity Vf, jvp is a decreasing function of 
the amplitude of the applied potential. In Fig. 2b, we 
plot the dependence of the inverse effective mass of the 



lower-band of the semi-Dirac point at the BZB (in the 
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vicinity of (k x , k y ) = ( — , 0)) on the applied potential pe- 

riod L and amplitude Uq respectively. We find that the 
effective mass along the y-direction changes sign from 
positive to negative as a function of L and Uq. The dis- 
persion at small k y for parameters U, Lq corresponding 
to negative effective mass is shown in Fig. 3(a). In this 
case, the bands at larger k y turn up leading to the emer- 
gence of another fermi surface. There is no such fermi 
surface doubling for parameters with positive mass. 



III. OTHER FORMS OF SPIN DENSITY WAVE 



A. The HSDW with U z / U y 

In this subsection, we shall discuss the case of a HSDW 
with different amplitudes in the y and z directions (i.e. 
U y 7^ V z) with the propagation direction of the spin den- 
sity wave still along the i direction. In order to study the 
effects of helical spin density wave with U z ^ U y , we first 
use the perturbation theory to calculate the renormalized 
velocity in the BZ center. In the limit of a perturbatively 
weak coupling to the HSDW < 1), the renormal- 

ization of group velocity (near the Dirac cone) given in 
Eq. [6j within second order perturbation theory should be 
modified to: 



Vf 



V- 



87T 2 h 2 V 2 F 



[2C/2 + [/f(l + cos20 k )] (11) 



It is clear that the above formula goes to Eq. [6jfor the 
proper HSDW. We notice that the spin density wave in 
the y direction does not induce the anisotropic group 
velocity near the Dirac cone, which is also confirmed by 
the full numerical calculation. We present the angular 
dependence of the renormalized group velocity for two 
extreme limits U y = and U z = in Fig. [5] 

We shall next explain the conditions for the appear- 
ance of the Dirac cones on the BZ boundary, which are 
supported by full numerical calculations as presented in 
Fig. |B] and [7J The largest contribution to the energy 
eigenvalues at the BZ boundary comes from these two 
degenerate states |s,k) and |s,k— (2n/L,0)), which are 
degenerate before applying the periodic potential. Only 
considering the above two states and using the relation 
$k-(27r/L,o) = 7r — on the BZB, we can write down the 
two band Hamiltonian within perturbation theory as: 



H(k) 



I Cjf + k 2 



\{sU y — U z cos0k)(«cos6'k + sin#k) 



V 



\{sU y — U z cos6>k)(— i cos6>k + sin^k) 



shv F Mj) 2 + k 2 



(12) 
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where s = ±1 denote the band index and 6*k is the angle 
between the vector k and the x direction. Diagonaliz- 
ing the above Hamiltonian we could get two eigenvalues 
written as: 




AE =| sU y - U z cos6 k |, 

From the last equation in Eq. I13[ we can clearly see that 
there are gapless states when sll y = U z cos#k, since < 
cos 0k < 1 for k on the BZ boundary (i.e. k = (tv/L, k y )). 
Therefore, we can have gapless states only for U z > U y . 
In addition, the energy charge carrier dispersion becomes 
parabolic along both x and y directions, which can be 
seen from Fig. EJa) and[7ta). For the clockwise helix, the 
energy gap opened on the BZ boundary near k = (w/L,0) 
is shown as a function of k y for three different values 
of U zy . In Fig. IDJb) and Fig. [TJb) the energy gap is 
plotted as a function of k y for U y /U z = 0.8 and U z /U y = 

0. 8 respectively. It is seen that the energy gap is non- 
monotonic function of \k y \ in the former, and monotonic 
in the latter case. The helical spin density with U y ^ 
U z affect the energy charge carrier dispersion on the BZ 
center similar to the proper HSDW, rendering the Dirac 
cone on the BZ center highly anisotropic. The velocity 
of the Dirac cone along x direction (the helix propagates 
along the x direction) is smaller than that perpendicular 
to x. When we go to the extreme limit , where one of the 
components U VtZ is set to zero (the spin density wave), 
the energy gap in the presence of the spin density wave 
with magnetization in only one direction {y or z), further 
simplifies to: 

f | E/ 2 cos0 k I, u y = o 
AE = \ (14) 
l\U y \, U z = 

Full numerical calculation shows (see Fig [5]) the depen- 
dence on k y for the energy gap opened on the BZB with 
the presence of the spin density waves in either z or y di- 
rection. The latter dependence is not predicted correctly 
within the second order perturbation theory which sug- 
gests that the energy gap for U z — is independent of 
k y (see second line of Eq. [T4")) . Finally, one should note 
that the Hamiltonian has different symmetries along y 
and z directions. For a spin density wave only in the 
y direction, the total Hamiltonian satisfies the relation: 
a z H<j z = —H, from which we can see that the upper 
and the lower band of the Dirac cone arc symmetric. 
While for the spin density wave only in the z direction, 
the total Hamiltonian satisfies particle-hole symmetry, 

1. e. a x H*(— h)a x = —H(k). The upper and the lower 
band of the Dirac cone are also symmetric when the spin 
density wave only has component in z direction. The ve- 
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FIG. 5: (Color online) (a) The group velocity Vf, (measured 
from the Dirac point) of charge carriers on the surface of the 
TI in units of the Fermi velocity {vf) versus the angle (9k 
with L = 90nm and U y = 0. Dotdashed red, solid green and 
dashed blue lines correspond to U z being 10 meV, 30 meV and 
50 meV, respectively, (b) The group velocity Vf, (measured 
from the Dirac point) of charge carriers on the surface of the 
TI in units of the Fermi velocity (vf) versus the angle (9k 
with L = 90nm and U z = 0. Dotdashed red, solid green and 
dashed blue lines correspond to U y being 10 meV, 30 meV 
and 50 meV, respectively. 
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FIG. 6: (Color online) (a) Energy charge carrier dispersion 
versus the wavevector k with k measured from (k x ,k y ) — 

TV 

(y,0) on the BZB. Dashed (green) and solid (blue) lines 

show the linear dispersion along the x direction and the 
quadratic dispersion along the y direction, respectively with 
U z = 50meV, U y = 40meV, L = 90nm and Go = 27r/90nm" 1 
(E(tv/L,0) was shifted to zero). (b)The energy gap AE be- 
tween the first and the second band at the BZB versus k y 
for charge carriers above the Dirac point with L = 90nm and 
U y — 0.8U Z . Dotdashed (red) , solid(green) and dashed (blue) 
lines correspond to U z being 10 meV, 30 meV and 50 meV, 
respectively. 



locity perpendicular to the x, which is the same as the 
original Fermi velocity of B^Ses, remains unaffected in 
the presence of the periodic potential in the z direction. 



B. Cycloidal spin density wave 

In this subsection, we analyze the effects of the cy- 
cloidal spin density wave on the surface states of topo- 
logical insulator. In the presence of the cycloidal spin 
density wave on top of the topological insulato r 31 ' 32 ; the 
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FIG. 7: (Color online) (a) Energy charge carrier dispersion 
versus the wavevector k with k measured from (k x ,k y ) = 

7T 

(y,0) on the BZB. Dashed (green) and solid (blue) lines 

Li 

show the linear dispersion along the x direction and the 
quadratic dispersion along the y direction, respectively with 
U z = 40meV, U y = 50meV, L = 90nm and Go = 27r/90nm" 1 
(E(n /L,0) was shifted to zero). (b)The energy gap AE be- 
tween the first and the second band at the BZB versus k y 
for charge carriers above the Dirac point with L — 90nm and 
U z = 0.8U y . Dotdashed (red) , solid(green) and dashed (blue) 
lines correspond to U y being 10 meV, 30 meV and 50 meV, 
respectively. 



a b 
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FIG. 8: (Color online) (a) The energy gap AE between the 
first and the second band at the BZB versus k y for charge 
carriers above the Dirac point with L = 90nm and U y = 0. 
Dotdashed (red) , solid (green) and dashed (blue) lines cor- 
respond to U z being 10 meV, 30 meV and 50 meV, respec- 
tively. (b)The energy gap AE between the first and the sec- 
ond band at the BZB versus k y for charge carriers above the 
Dirac point with L — 90nm and U z = 0. Dotdashed (red) , 
solid(green) and dashed (blue) lines correspond to U y being 
fO meV, 30 meV and 50 meV, respectively. 



potential is modified to: 

2n 2n 
U c (x) = U x a x cos(— x) + Uy. z <Ty. z sin(— x), (15) 

Li Li 

The term involving a x in the magnetic potential will dis- 
appear under the gauge transformation: 

2-7T 

e -iA( x ) He iA{ X ) = Hq + Uy z(Jy z sin ( — x), (16) 

where Hq is the single Dirac fermion Hamiltonian with- 
out applied potential as given in Eq. [T] and the function 



A(x) is given by: 

A(x) = — — sm(— ) 17 

For the general periodic function of spin density wave 
along x direction, i.e. F(x)a x , the a x term will disap- 
pear if we apply the gauge transformation e lA p( x ) with 
Ap(x) = — J F(x')/hdx' (c is any constant). Thus, the 
effects of the cycloidal spin density wave comes from the 
term with <Jy lZ , which goes to the extreme limit as dis- 
cussed in Sec lHI Al 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have considered the effects of a helical 
spin density wave on the surface states of a topological 
insulator and also the effects of other spin density wave 
on the surface transport property of topological insula- 
tor. We find that the HSDW acts like an effective spin 
potential on the TI surface and breaks the local time 
reversal invariance in the latter. The applied spin poten- 
tial has two main consequences: First, the group velocity 
of electrons at the Dirac point is strongly suppressed in 
the direction transverse to the periodic potential. Sec- 
ondly, new semi-Dirac points emerge for the first upper 
(lower) band at the Brillouin zone boundaries that corre- 
sponds to the right-handed (left-handed) chirality of the 
applied proper HSDW. The semi-Dirac points are char- 
acterized by linear dispersion parallel to the direction x 
of the applied periodic potential but quadratic dispersion 
perpendicular to the x direction. When the chemical po- 
tential of the TI is such that the semi-Dirac points at the 
BZBs give the main contribution to transport, wc expect 
the surface transport properties to be highly anisotropic. 
This is because at these points, the dispersion is linear in 
one direction and quadratic in the other. The momentum 
space location of these new semi-Dirac cones can also be 
manipulated by applying an in-plane magnetic field to 
the HSDW. Such a magnetic fielcU^ can change the ori- 
entation of the spin rotation axis relative to the pitch 
vector of the HSDW. By studying the effects of such a 
change in the applied periodic potential on the new semi- 
Dirac cones and measuring the anisotropic component of 
the transport, it should be possible to isolate the surface 
state contribution from the bulk in the total conductance 
in topological insulators with relatively small bulk gaps. 
The number of Dirac cones appearing on the BZ bound- 
ary is determined by the relative ratio of U z /U y , which 
could be zero, one or two corresponding to U z /U y < 1, 
Uz/Uy = 1 and U z /U y > 1, respectively. We also prove 
that the spin density wave with spin vectors along x di- 
rection will not change the surface transport property 
of a three dimensional topological insulator. One par- 
ticular advantage of our proposal is that the proposed 
experiment can be carried out in the existing TI systems 
without any problem arising from the bulk conduction (as 
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a result of unintentional bulk doping by defects), which is 
invariably present in almost all current TI systems, since 
the bulk transport is isotropic and presumably unaffected 
by the presence of the HSDW. 

Although the proposed system in this work is a sand- 
wich structure of a TI and a multiferroic, we believe that 
this structure may turn out to be a suitable candidate 
for the direct manifestation of the elusive 2D transport 
properties even in the presence of considerable bulk con- 
duction since the anisotropy introduced by the presence 
of the HSDW would only affect the surface transport 
properties without affecting much the bulk conduction 
behavior. Given the great current interest and activity 
in the observation of the 2D surface transport properties 
in TI materials^—, we are optimistic that our proposed 
structure could go a long way in establishing the exper- 



imental behavior of 2D surface transport in topological 
insulators. 
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